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Abstract 

In  this  paper,  a  new  strategy  for  scaling  burners  based  on  “mild  combustion”  is  evolved  and  adopted  to 
scaling  a  burner  from  3  to  a  150  kW  burner  at  a  high  heat  release  rate  of  5  MW/m3.  Existing  scaling  meth¬ 
ods  (constant  velocity,  constant  residence  time,  and  Cole’s  procedure  [Proc.  Combust.  Inst.  28  (2000) 
1297])  are  found  to  be  inadequate  for  mild  combustion  burners.  Constant  velocity  approach  leads  to 
reduced  heat  release  rates  at  large  sizes  and  constant  residence  time  approach  in  unacceptable  levels  of 
pressure  drop  across  the  system.  To  achieve  mild  combustion  at  high  heat  release  rates  at  all  scales,  a  mod¬ 
ified  approach  with  high  recirculation  is  adopted  in  the  present  studies.  Major  geometrical  dimensions  are 

t  /o  _ 

scaled  as  D  ~  Q  with  an  air  injection  velocity  of  ~100  m/s  (Ap  ~  600  mm  water  gauge).  Using  CFD  sup¬ 
port,  the  position  of  air  injection  holes  is  selected  to  enhance  the  recirculation  rates.  The  precise  role  of 
secondary  air  is  to  increase  the  recirculation  rates  and  burn  up  the  residual  CO  in  the  downstream.  Mea¬ 
surements  of  temperature  and  oxidizer  concentrations  inside  3  kW,  150  kW  burner  and  a  jet  flame  are  used 
to  distinguish  the  combustion  process  in  these  burners.  The  burner  can  be  used  for  a  wide  range  of  fuels 
from  LPG  to  producer  gas  as  extremes.  Up  to  8  dB  of  noise  level  reduction  is  observed  in  comparison  to 
the  conventional  combustion  mode.  Exhaust  NO  emissions  below  26  and  3  ppm  and  temperatures  1710 
and  1520  K  were  measured  for  LPG  and  producer  gas  when  the  burner  is  operated  at  stoichiometry. 
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1.  Introduction 

Guidelines  for  scaling  are  important  in  the  de¬ 
sign  of  combustion  systems.  A  large  number  of 
dimensionless  groups  for  scaling  are  proposed  by 
Spalding  [1]  and  Beer  and  Chigier  [2].  It  is  recog¬ 
nized  that  maintaining  all  the  variables  constant 
during  the  process  is  not  possible  partly  because 
of  internal  inconsistencies.  It  is  imperative  to 
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adapt  scaling  based  on  a  selected  set  of  non-di¬ 
mensional  quantities. 

Constant  velocity,  CV  and  constant  residence 
time,  CRT  approaches  have  been  applied  to  scale 
up  the  burners  and  furnaces  from  laboratory  scale 
[3-10].  For  a  burner,  the  total  thermal  input  is  gi¬ 
ven  as  Q  =  rh{H  =  KpU0D20.  In  CV  approach,  the 
burner  inlet  velocity  is  maintained  constant,  and 
geometrical  dimensions  are  derived  from  the  rela¬ 
tionship  DfiDi  =  (QilQi)112-  For  CRT  approach, 
the  ratio  D0/U0  (inertial  or  convective  timescale) 
is  maintained  constant  while  increasing  the  burner 
thermal  input.  The  new  physical  dimensions  are 
determined  through  a  relationship  D2!D\  -  (Q2l 


1540-7489/$  -  see  front  matter  ©  2004  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved, 
doi:  10.1016/j  .proci.2004.07 .045 


2614 


S.  Kumar  et  al.  I  Proceedings  of  the  Combustion  Institute  30  ( 2005 )  2613-2621 


Nomenclature 

CRT 

Constant  residence  time  scaling 

Q,  Qu  Q2 

Burner  thermal  input  (kW) 

approach 

U0,  ua,  u{ 

Inlet  velocity  of  air/fuel  or 

CV 

Constant  velocity  scaling  approach 

characteristic  velocity  (m/s) 

D,D0, 

Different  burner  dimensions 

tha,  fhf 

Air  and  fuel  flow  rates 

Di,D2 

Q 

Heat  release  rate  per  unit  volume 

/v 

Volume  fraction 

MW  /m3 

H 

Calorific  value  of  fuel 

P 

Density  of  air/fuel  (kg/m3) 

LPG 

Liquefied  petroleum  gas  (~80% 

butane  and  20%  propane) 

v  IIllAlllg} 

%  ?  % 

Characteristic  mixing  time  or 

Producer 

gas 

Low  calorific  value  gas  (LL 
-  20%,  CO  -  20%,  CHa  -  2%, 

CO2  ~  13%  and  rest  N2) 

convective  timescales,  D0/U0  (Is) 

i  /'■y 

Q i)  .  In  another  approach  adopted  by  Cole  et  al. 
[8],  both  air  velocity  and  jet  area  are  increased 
equally  for  scaling  combustors.  The  velocity  and 
jet  diameters  are  scaled  as  U2IUi  =  (Q2/Q\)1  2  and 
D2IDi  -  (Q2IQ\)114  to  burner  inputs.  More  details 
about  the  variation  of  critical  parameters  with  dif¬ 
ferent  scaling  approaches  are  given  in  Table  1. 

It  is  known  that  CV  approach  increases  the 
characteristic  mixing  time  and  reduces  the  rate 
of  mixing  [3-10].  To  maintain  constant  rate  of 
mixing,  the  inlet  velocity  should  be  increased  as 
Ql/ 3  with  burner  thermal  input  [4-6,8].  In  Cole’s 
[8]  approach,  jet  velocity  increases  at  faster  rate 
than  jet  diameter.  Therefore,  both  CRT  and 
Cole’s  [8]  approaches  lead  to  large  pressure  drop 
across  the  combustion  system.  The  experimental 
investigations  using  Cole’s  [8]  approach  on  an 
acoustically  excited  combustor  showed  consistent 
performance  for  pollutant  emissions,  flame  stabil¬ 
ity,  and  enhanced  mixing  at  smaller  levels.  The 
improvement  in  emissions’  performance  at  larger 
power  scales  is  reported  to  be  insignificant. 


The  scaling  studies  on  swirl  stabilized  pulverized 
coal  burners  have  shown  that  NOx  emissions  de¬ 
pend  on  local  fluid  flow  behavior  in  the  internal 
recirculation  zone  [4,5].  Computational  investiga¬ 
tions  by  Weber  and  Breussin  [6]  predicted  that  be¬ 
yond  a  certain  thermal  input,  NO*  emissions 
remained  independent  of  the  scaling  approach  used. 

CV  approach  fails  to  produce  aerodynamic  sim¬ 
ilarity  in  the  near  burner  region  for  swirl  stabilized 
natural  gas  burners  [7].  This  is  a  critical  factor  for 
NOx  formation  in  the  gas  burners.  Detailed  analy¬ 
sis  of  NOt  emissions  from  two  gas  burners  at  67  and 
266  kW  thermal  levels  showed  the  importance  of 
prompt  NO  formation  in  the  near  burner  zone  of 
a  combustion  system  [10]. 

The  exhaust  gas  recirculation  for  NOx  reduc¬ 
tion  from  combustion  systems  has  drawn  interest 
due  to  its  promising  features  [11-16].  When  recir¬ 
culation  rates  are  high  enough  and  at  tempera¬ 
tures  greater  than  the  auto-ignition  temperature 
of  the  fuel,  a  stable  combustion  mode  exists.  This 
combustion  mode  is  known  as  mild  or  flameless 


Table  1 


Comparison  of  various  parameters  in  different  scaling  approaches 


Scaling  approaches 

Geometric  scaling  D  =  (D2ID{) 

Velocity  scaling  U  =  {U2IU\) 

r  mixing 

Re 

q'" 

CV 

~Q1/2 

Constant 

~Q1/2 

~Q1/2 

~Q-1/2 

CRT 

~Q1/3 

~QI/3 

Constant 

~q2/3 

Constant 

Cole  [9] 

~Q1/4 

~QI/2 

~Q-1/4 

~q3/4 

~QI/4 

Present 

~Q1/3 

~100  m/s 

~Q1/3 

~Q1/3 

Constant 

Table  2 

Summary  of  the  previous  work  in  mild  combustion  and  residence  times  used  in  these  experiments 

Ref. 

Uf  (m/s) 

Tf  (ps) 

Ca  (m/s) 

Ta  (|is) 

q'"  ( MWIrrii ) 

Q  (kW) 

[12] 

20 

250 

73.7 

74.6 

0.32 

10 

[13] 

9.34 

503.2 

33 

151.51 

0.18 

6 

[14] 

12.57 

318.21 

28.9 

162.58 

0.18 

6 

[18] 

100 

100 

70 

1771 

0.023 

580 

[20] 

7.9-70.7 

114-4.2 

— 

— 

— 

— 

[15] 

20-100 

25-5 

26-130 

77-15.5 

5.6 

1-5 

Present 

243 

3 

95 

52 

5.6 

150 
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combustion  [11-19].  Table  2  summarizes  the  work 

carried  out  on  mild  combustion.  Lower  convective 

timescales  seem  to  be  a  critical  factor  to  achieve 

•  ///  2 

mild  combustion  with  high  Q  (  5  MW/m  )  and 
reactants  at  ambient  temperature.  The  previous 
experiments  are  conducted  in  the  thermal  range 
of  1-580  kW  with  low  heat  release  rates  (23- 
320  kW/m3)  [12-15]. 

To  scale  a  burner  with  high  Q  ,  if  one  uses  a 
CRT  method,  jet  velocity  increases  as  Q  and 
hence  leads  to  unacceptable  levels  of  pressure 
drop  across  the  system  beyond  a  certain  thermal 
input  range  as  shown  in  Table  3.  One  needs  to  ex¬ 
plore  other  alternatives  to  scale  a  mild  combus¬ 
tion  burner  while  maintaining  the  geometric, 
dynamic,  and  thermal  similarities. 

The  objectives  of  the  current  research  are  to 
use  the  results  of  3  kW  laboratory  scale  burner, 
scale  it  to  a  large  level,  in  this  case  150  kW,  estab¬ 
lish  mild  combustion  in  a  high  heat  release  burner, 
and  suggest  scaling  laws  for  the  mild  combustors. 

2.  Computations 

The  objectives  of  the  computational  studies  are 
to  optimize  the  burner  geometry,  to  quantify  the 
recirculation  rates  and  to  predict  the  combustion 
and  fluid  flow  behavior  of  a  150  kW  mild  combus¬ 
tion  burner.  The  same  code  that  was  used  for 
3  kW  laboratory  burner  is  used  here.  The  details 
related  to  computational  strategy,  fluid  flow,  and 
combustion  modeling  of  the  burner  can  be  found 
in  Sudarshan  et  al.  [15].  Since  the  geometry  pre¬ 
sents  a  sixfold  symmetry  with  six  alternate  fuel 
and  air  injection  jet  arrangements  along  the  cen¬ 
tral  axis,  one-sixth  part  of  the  burner  is  considered 
for  the  numerical  simulation.  To  obtain  grid  inde¬ 
pendent  results,  grid  resolution  studies  are  carried 
out  with  the  number  of  grid  points  varying  from 
100,000  to  1,000,000.  The  results  with  respect  to 
600,000  grid  points  were  within  1%  for  all  meshes 
up  to  1,000,000  grid  points. 

3.  Geometry  optimization  and  computational  results 

The  3  kW  laboratory  burner  investigated  earlier 
by  Sudarshan  et  al.  [15]  is  scaled  using  CV,  CRT, 


and  Cole’s  [8]  scaling  principles.  Table  3  shows 
the  detailed  dimensions,  velocities,  and  other  re¬ 
lated  details  of  the  1 50  kW  scaled  burner  with  dif¬ 
ferent  approaches.  The  burner  is  theoretically 
scaled  to  2  MW  to  show  the  effect  of  different  ap¬ 
proaches  on  recirculation  rates  and  heat  release 
rates.  Figure  1  shows  the  effect  of  different  scaling 
approaches  on  physical  geometry  of  150  kW  bur¬ 
ner.  The  physical  dimensions  of  the  scaled  burner 
decrease  from  CV  to  CRT  and  Cole’s  [8]  approach, 
and  corresponding  Q  increases.  Figure  2  shows 
the  variation  of  recirculation  rates  in  1 50  kW  scaled 
burner  with  different  scaling  approaches  obtained 
from  computational  studies.  Curve  (a)  shows  that 


Air  and  fuel  inlets 


Fig.  1.  Details  of  the  150  kW  scaled  burner  using 
different  approaches. 


Table  3 


Summary  of  the  geometrical  dimensions  with  different  scaling  approaches 


Scale  factor 

1 

10 

50 

667 

Q  (kW) 

3 

30 

150 

2000 

Z)a,  £/a  for  CRT 

2,  80 

4.3,  172 

7.4,  291 

17.5,  698 

Dx,  D2,  L  (mm)  for  CRT 

60,  90,  120 

130,  194,  258 

221,  331,  440 

525,  790,  1050 

Di,  D2,  L  (mm)  for  CV 

60,  90,  120 

190,  285  380 

425,  636,  850 

1550,  2325,  3100 

D\,  D2,  L  (mm)  Cole  [8] 

60,  90,  120 

107,  160,  214 

160,  240,  320 

305,  458,  610 

Present  Da,  Ua 

2,  80 

— 

5,  95 

— 

Present  ia  (|is) 

25 

— 

52.4 

— 

Recirculation  rate 

2.8 

— 

2.3 

— 
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Thermal  Input  (kW) 


Fig.  2.  Variation  of  recirculation  rate  and  heat  release 
rates  with  thermal  power.  Curve  (a)  Recirculation  rate 
variation  with  CRT  approach.  Curve  (b)  Recirculation 
rate  variation  with  CV  approach.  Curve  (c)  Recircula¬ 
tion  rate  variation  when  major  burner  dimensions  are 
determined  through  D^Q113  with  an  air  injection  veloc¬ 
ity  of  100  m/s.  Curve  (d)  Variation  of  heat  release  rate 
for  CRT  approach.  Curve  (e)  Variation  of  heat  release 
rate  for  CV  approach. 


the  recirculation  rates  drop  from  280%  to  220%  as 
the  burner  is  scaled  from  3  kW  level  to  2  MW  using 
CRT  approach.  The  air  inlet  velocity  (scaled 
as  UQ  ~  Q1/3)  increases  from  79  to  698  m/s.  This 
leads  to  unacceptable  level  of  pressure  drop  across 
the  combustion  systems.  Curve  (d)  shows  that  cor¬ 
responding  heat  release  rate  remains  constant  with 
CRT  approach.  Similarly  Cole’s  [8]  (U0~Ql/2)  ap¬ 
proach  also  results  in  large  pressure  drop  across 
the  combustion  systems. 

When  CV  scaling  approach  is  used  to  deter¬ 
mine  the  burner  dimensions,  the  recirculation 
rates  drop  from  280%  to  190%  as  shown  by  curve 
(b).  Curve  (e)  represents  the  corresponding  Q 
variation  with  CV  approach.  The  heat  release 
rates  drastically  drop  as  1  IQ112  from  5.6  to 
0.217  MW/m3  as  burner  is  scaled  from  3  kW  to 
2  MW.  Curve  (c)  represents  variation  of  recircula¬ 
tion  rates  in  the  burner  when  the  major  burner 

.  .  •  /// 

dimensions  are  determined  to  maintain  high  Q 
at  5.6  MW/m3  and  air  inlet  velocity  at  ~100  m/s, 
an  affordable  choice  in  industrial  applications. 
To  maintain  Q  constant,  combustion  system  vol¬ 


ume  should  be  increased  in  proportion  to  thermal 
input.  Therefore,  major  burner  dimensions  are 
scaled  as  D  ~  Ql/3.  In  the  current  scaling  strategy, 
recirculation  rates  drop  to  153%  at  150kW  and 
70%  at  2  MW  level,  thus  making  it  difficult  to 
achieve  mild  combustion  in  the  scaled  burners. 

At  this  point  of  time,  a  150  kW  burner  scaled 
from  3  kW  laboratory  burner  is  tested  experimen¬ 
tally  for  the  demonstration  of  mild  combustion 
mode.  The  major  dimensions  are  determined 
using  D  ~  Q  to  maintain  high  Q  .  Air  jet  veloc¬ 
ities  of  100  m/s  are  considered.  This  experimental 
burner  is  tested  with  both  LPG  and  producer  gas 
(typical  producer  gas  composition  CO  20%,  H2 
20%,  C02  13%,  CH4  2%,  and  rest  N2).  It  is  ob¬ 
served  that  due  to  low  recirculation  rates  and 
large  air  and  fuel  jet  diameters  (large  convective 
timescales  of  air  and  fuel  jets),  the  combustion 
zones  are  clearly  visible  as  a  kind  of  highly  con¬ 
fined  jet  flames  attached  to  either  air  or  fuel  jets 
(for  detailed  operating  conditions  see  Table  4). 
Overall  observed  emission  levels  are  low.  The 
presence  of  highly  confined  jet  flames  in  the  com¬ 
bustion  zone  prompted  further  investigations. 

The  steep  reduction  in  recirculation  rates  led  to 
the  exploration  of  other  alternatives  to  increase 
the  overall  recirculation  rates  to  achieve  mild 
combustion.  Initial  trials  with  different  injection 
schemes  showed  that  various  air/fuel  injection 
schemes  had  very  little  effect  on  overall  recircula¬ 
tion  rates.  The  air  and  fuel  injection  details  are 
shown  in  Fig.  3.  Both  air  and  fuel  are  injected 
as  a  set  of  six  holes  at  different  locations  at 
90  mm  pitch  diameter.  The  diameters  and  number 
of  air  and  fuel  jets  are  estimated  to  keep  the  DJU0 
ratio  in  the  same  range  as  for  laboratory  scale 
burner  operational  range  (see  Table  4).  The  recir¬ 
culation  rates  are  further  enhanced  by  appropri¬ 
ately  optimizing  the  position  of  secondary  air 
injection  (~20%  of  total  air).  The  secondary  air 
is  injected  as  a  set  of  multiple  high-speed  jets  from 
the  top.  A  number  of  calculations  are  carried  out 
to  reveal  the  effect  of  location  and  velocity  of  sec¬ 
ondary  air  injection  on  recirculation  rates.  Figure 
4  shows  the  variation  of  the  recirculation  rate  with 
the  position  of  secondary  air  with  respect  to  the 


Table  4 


Summary  of  experiments  carried  out  on  the  mild  combustion  burners  (MC-mild  combustion,  AF-attached  flames,  and 
PG-Producer  gas) 


Burner 

Q  (kW) 

Uf  (m/s) 

Df  (mm) 

Tf  (|IS) 

Fa  (m/s) 

Da  (mm) 

Ta  (|1S) 

Fuel 

Remarks 

I 

1 

20 

0.5 

25 

27 

2 

74 

LPG 

MC 

3 

60 

0.5 

8.3 

80 

2 

25 

LPG 

MC 

5 

100 

0.5 

5 

135 

2 

14.8 

LPG 

MC 

II 

150 

125 

2 

16 

100 

10 

100 

LPG 

AF 

150 

95 

12 

126 

128 

10 

78 

PG 

AF 

150 

63 

6 

95 

78 

5 

64 

PG 

AF 

150 

92 

5 

54 

78 

5 

64 

PG 

MC 

150 

243 

0.7 

2.9 

95 

5 

52.4 

LPG 

MC 
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Fig.  3.  Details  of  optimized  configuration  for  150  kW 
burner  with  alternate  peripheral  injection  schemes  for 
both  LPG  and  producer  gas  fuels. 


and  02-volume  behavior  is  extracted  from  the  cal¬ 
culations.  The  volume  elements  for  AT  =  100  K 
and  AX0,  =  0.01  steps  are  determined  over  the  en¬ 
tire  combustor  and  plotted  for  (a)  a  turbulent  jet 
diffusion  flames  (b)  1 50  kW  optimized  burner, 
and  (c)  3  kW  laboratory  burner.  Figures  5  and  6 
show  the  cumulative  distribution  of  volume  frac¬ 
tion  variation  with  temperature  and  02  for  burners 
at  different  scales.  For  a  turbulent  jet  flame, 
/v  =  0.54  for  temperature  <1000  K  and  /v  =  0.54 
for  02  mass  fraction  >0.15.  For  150  kW  optimized 
burner,  /v  =  0.93  for  temperature  >1000  K  (auto¬ 
ignition  temperature  of  the  fuel)  and  is  almost  uni¬ 
formly  distributed  over  the  whole  range.  Similarly, 
/v  =  0.96  for  02  mass  fraction  <0.15  and  is  uni¬ 
formly  distributed  in  the  range  of  0-0.15.  For  a 
3  kW  mild  combustion  burner,  fy  =  0.93  for  tem¬ 
perature  >1300  K  and  for  02  mass  fraction  <0.07, 
/v  =  0.97.  This  deviation  of  temperature  and  oxi¬ 
dizer  mass  fraction  distributions  from  3  kW  mild 
combustion  burner  [15]  appears  significant.  This 
is  attributed  to  the  fact  that  in  the  current  burner, 
the  number  of  air-fuel  jets  is  six  times  larger  than 
the  previously  investigated  3  kW  burner.  The  group 


Secondary  air  injection  velocity  (m/s) 

0  20  40  60  80  100  120 


Fig.  4.  Variation  of  recirculation  rate  with  secondary  air 
position  and  velocity.  Curve  (A)  Variation  of  recircula¬ 
tion  rate  with  secondary  air  injection  position.  Curve  (B) 
Variation  of  recirculation  rate  with  secondary  air 
velocity  at  optimum  position. 


wall  and  secondary  air  velocity.  Curve  (a)  shows 
that  the  position  of  secondary  jets  has  a  strong  ef¬ 
fect  on  the  recirculation  rate.  The  recirculation 
rate  reaches  a  maximum  of  196%  for  a  constant 
injection  velocity  of  24  m/s  when  injection  holes 
are  located  at  1 1  mm  from  the  wall.  Curve  (b) 
shows  that  recirculation  rate  varies  almost  linearly 
with  the  secondary  injection  velocity.  Compared 
to  the  conventionally  injected  secondary  air  [15] 
at  1 50  kW  level,  the  recirculation  rates  are  en¬ 
hanced  by  80%. 

To  establish  a  quantitative  comparison  between 
the  burners  at  different  scales,  temperature-volume 


Fig.  5.  Predicted  cumulative  temperature-volume 
behavior  for  (a)  classical  turbulent  jet  flame,  (b) 
150  kW  scaled  burner,  (c)  3  kW  laboratory  scale  burner. 


O2  mass  fraction 


Fig.  6.  Predicted  cumulative  02-volume  behavior  for  (a) 
classical  turbulent  jet  flame,  (b)  1 50  kW  scaled  burner, 
(c)  3  kW  laboratory  scale  burner. 
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of  air  jets  influenced  a  large  volume  when  compared 
to  a  single  air  jet.  For  the  mild  combustion  mode, 
02  is  typically  in  the  range  of  0-0. 1 5,  and  tempera¬ 
ture  is  greater  than  1000  K  (auto-ignition  tempera¬ 
ture  of  the  fuel)  [17].  Hence,  even  case  (b)  can  be 
considered  as  in  mild  combustion  mode. 


4.  Scaling  of  mild  combustion  burners 

Table  4  shows  the  summary  of  convective  time- 
scales  for  different  combinations  of  air  and  fuel 
jets  employed  during  the  experimental  investiga¬ 
tions.  Highly  confined,  attached,  and  fluctuating 
jet  flames  appeared  in  the  combustion  zone  below 
1  kW.  The  reduction  in  the  air  and  fuel  flow  rates 
leads  to  reduction  in  recirculation  rates  and  re¬ 
sults  in  the  appearance  of  attached  flames  within 
the  reaction  zone.  At  this  operating  condition, 
convective  timescales  for  air  and  fuel  jets  are  74 
and  25  ps,  respectively. 

At  150  kW  thermal  level,  experiments  are  car¬ 
ried  out  with  different  air-fuel  injection  combina¬ 
tions  for  both  LPG  and  producer  gas.  Highly 
confined  jet  flames  are  observed  visibly,  attached 
to  either  air  or  fuel  jets  for  the  cases  of  convective 
timescales  greater  than  ~80  ps.  It  is  observed  that 
mild  combustion  is  obtained  successfully  with 
both  LPG  and  producer  gas  for  DJU0  ratio  below 
80  ps.  From  a  series  of  experiments  in  the  1- 
150  kW  range  on  mild  combustion  burners,  it  is 
concluded  that  for  successful  scaling  of  mild  com¬ 
bustion  burners  with  high  heat  release  rates,  DJ 
U0  ratio  should  be  maintained  below  ~80  ps. 

Mild  combustion  is  achieved  when  flames  are 
lifted  off  from  the  primary  burner  zone  at  high 
velocities.  This  can  be  explained  on  the  basis  of 
lift-off  concept  of  simple  jet  diffusion  flames, 
which  depends  on  local  temperature,  reactant 
concentration,  velocity,  and  diameter  of  the  injec¬ 
tion  jets.  Confined  jet  flames  are  expected  to  ap¬ 
pear  in  the  combustion  zone  as  long  the 
velocities  are  below  blow-off  and  recirculation 
rates  are  low. 


5.  Experiments 

Two  fuels  are  selected  for  the  experimental 
studies  on  150  kW  burner  to  show  that  the  burner 
can  be  used  over  a  wide  range  of  fuels.  Producer 
gas  and  LPG  are  chosen  as  extremes  (Calorific  va¬ 
lue  variation  4.5-45  MJ/kg).  The  burner  is  oper¬ 
ated  at  stoichiometry  with  1 50  kW  thermal 
input.  The  dimensions  of  the  burner  are  fixed  by 
the  total  thermal  input  and  Q  in  the  burner. 
These  dimensions  are  further  modified  through  a 
number  of  computations  aimed  at  optimizing  the 
burner  configuration.  The  details  of  the  air  and 
fuel  injection  schemes  for  LPG  and  producer  gas 
are  shown  in  Fig.  3.  Typical  air  and  LPG  mass 


flow  rates  are  50  and  3.2  g/s.  The  flow  rates  for 
producer  gas  and  air  are  39  and  47  g/s.  Eighty  per¬ 
centage  of  the  total  air  is  supplied  through  the  pri¬ 
mary  inlet  and  20%  through  the  secondary  inlets. 
The  temperatures  in  the  reaction  zone  are  mea¬ 
sured  by  using  50  pm  Pt-13%Pt-Rh  thermocou¬ 
ples.  Measured  temperatures  are  corrected  for 
heat  loss  by  radiation.  The  corrected  temperatures 
are  accurate  within  ±50  K  of  actual  temperature. 
Species  (NO,  CO,  C02,  and  02)  concentrations 
are  measured  in  the  reaction  zone  by  using  Quin- 
tox  KM-9106  flue  gas  analyzer.  A  specially  de¬ 
signed  water-cooled  stainless  steel  probe  is  used 
to  draw  the  sample  gases  from  the  reaction  zone. 
The  sample  gases  are  immediately  cooled,  dried, 
and  then  transferred  to  the  analyzer  continuously. 
A  Lutron  SL-4001  sound  level  meter  is  used  to 
measure  the  sound  levels  during  the  combustion 
experiments.  More  details  of  the  Quintox  gas  ana¬ 
lyzer  and  sound  level  meter  are  mentioned  in 
Sudarshan  et  al.  [15].  The  noise  level  measure¬ 
ments  are  taken  at  a  point  50  mm  from  outlet  of 
the  burner  at  the  same  plane. 


6.  Results  and  discussion 

Most  of  the  results  presented  in  this  section  are 
for  the  150  kW  case.  The  burner  is  operated  at  a 
stoichiometric  air/fuel  ratio.  Figure  7  shows  the 
distinction  between  the  conventional  combustion 
and  mild  combustion  operation  at  150  kW  level 
with  LPG  and  producer  gas.  The  conventional 
combustion  mode  is  achieved  by  reducing  the  air 
and  fuel  flow  rates  in  the  system  leading  to  lower 
recirculation  rates  and  shifts  in  operation  to  at¬ 
tached  flame  mode  as  in  Fig.  7 A.  The  burner  oper¬ 
ation  in  mild  combustion  mode  is  shown  in  Fig. 
7B-D.  The  injection  arrangement  is  clearly  visible 
and  totally  transparent  at  the  injection  plane.  A 
very  weak  flame  is  present  in  the  reaction  zone 
which  is  light  bluish  in  color  and  barely  visible. 
All  the  walls  are  red  hot  and  glowing  consistently. 
This  combustion  mode  is  achieved  by  large  recir¬ 
culation  rates  of  the  combustion  products  into 
the  fresh  reactants. 

The  acoustic  level  measurements  are  carried 
out  during  the  cold  flow  (only  air  jets),  mild  com¬ 
bustion  mode,  and  conventional  combustion 
mode.  The  measured  levels  are  103,  105,  and 
113  dB,  respectively.  Approximately  8  dB  of  noise 
reduction  is  observed  when  operational  mode 
shifts  from  conventional  combustion  to  mild  com¬ 
bustion  for  reasons  known  earlier  [15].  Similar 
noise  reduction  is  observed  when  the  burner  is 
operated  with  producer  gas. 

The  composition  of  species  is  measured  at  two 
axial  locations,  150  and  400  mm  downstream 
from  the  injection  plane.  Figure  8  shows  the  CO, 
C02,  and  02  mole  fractions  and  temperatures 
measured  at  150  mm  axial  position  across  section 
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Fig.  7.  Comparison  between  conventional  and  mild 
combustion.  (A)  Conventional  turbulent  combustion 
with  low  recirculation  rates.  (B,C)  Mild  combustion 
mode  with  LPG  fuel.  (D)  Mild  combustion  mode  with 
producer  gas  fuel. 
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Fig.  8.  Species  concentration  and  temperature  measure¬ 
ments  with  LPG  fuel  at  150  mm  from  injection  plane. 

A-A  (Fig.  3).  O2  mole  fraction  drops  from  10%  on 
air  jet  side  to  very  low  values  of  1%  on  fuel  jet 
side.  02  is  fairly  well  distributed  with  small  gradi¬ 
ents  across  the  measuring  plane.  Measured  02 
mole  fraction  clearly  indicates  the  air  and  fuel 
jet  injection  sides  (across  section  A-A).  The  spe¬ 
cies  composition  variation  across  this  plane  is 
moderately  small.  Low  species  gradients,  high 
temperature,  and  low  concentration  of  reacting 
species  suggest  the  presence  of  a  slow  reaction 
over  a  large  area.  The  measured  temperature  var¬ 
ies  in  the  range  of  1200-1550  K  and  fairly  uniform 
over  the  reaction  zone. 

Figure  9  shows  the  species  and  temperature 
measurements  400  mm  downstream.  The  temper¬ 
atures  are  far  more  uniform  and  vary  between 
1500  and  1750  K.  The  temperature  gradients  at 
this  plane  are  much  smaller  than  those  compared 
to  150  mm.  The  species  concentration  variation  at 
this  plane  is  very  small.  The  uniformity  in  species 
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Fig.  9.  Species  compositions  and  temperature  measure¬ 
ments  with  LPG  fuel  at  400  mm  from  injection  plane. 

composition  across  the  plane  suggests  at  the  con¬ 
tinuance  of  slow  combustion  reaction  at  this 
plane.  At  exhaust,  the  emissions  recorded  are 
26  ppm  NO,  1%  CO,  and  an  average  temperature 
1710  K.  The  CO  emissions  are  in  the  range  of  pre¬ 
viously  reported  experiments  [15].  It  is  observed 
during  the  experiments  that  very  low  CO  emis¬ 
sions  (~0.03%)  are  recorded  when  ~10%  more 
air  is  added  downstream  to  dilute  the  combustion 
products  and  burn  the  residual  CO. 

The  injection  plane  of  the  burner  is  slightly 
modified  to  operate  the  same  burner  with  producer 
gas  (calorific  value  ~4.5  MJ/kg)  as  shown  in  Fig. 
3B.  The  typical  stoichiometric  ratio  for  producer 
gas  is  ~1.2.  A  200  kW  thermal  level  woody  bio¬ 
mass-based  gasifier  continuously  supplied  pro¬ 
ducer  gas  for  burner  operation  [21,22].  Figure  10 
shows  the  species  and  temperature  measurements 
at  400  mm  across  section  A-A.  The  temperatures 
measured  at  400  mm  are  quite  uniform  across  the 
radial  plane.  The  important  point  to  note  is  that 
temperature  is  much  more  uniform  in  the  reaction 
zone  than  compared  to  the  LPG  fuel.  Temperature 
variations  across  the  radial  plane  are  less  than 
200  K,  and  mean  temperature  is  above  1400  K. 
The  average  temperature  measured  across  the  ex¬ 
haust  plane  is  1 520  K.  Large  concentrations  of 
02  on  one  side  and  CO  on  the  other  side  indicate 
the  approximate  position  of  the  air  and  fuel  jets. 
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Fig.  10.  Species  concentration  and  temperature  mea¬ 
surements  with  producer  gas  fuel  at  400  mm  from 
injection  point. 
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Small  gradients  of  CO,  C02,  and  02  describe  the 
distributive  and  sluggish  nature  of  reaction  zone 
over  a  large  area.  Similar  behavior  of  species  con¬ 
centration  and  temperatures  are  recorded  at 
150  mm  from  the  injection  plane.  In  contrast  to 
the  LPG  operated  burner,  the  NO  emissions  from 
producer  gas  operation  are  very  low.  The  measured 
CO  and  NO  emissions  at  the  exhaust  are  0.211% 
and  3  ppm  against  the  1%  and  26  ppm  for  LPG. 
The  difference  in  emissions  could  be  attributed  to 
a  difference  in  the  average  operating  temperature 
(~200  K)  and  in  calorific  values  of  two  fuels. 

7.  Summary 

The  proposed  scaling  approach  is  shown  to  be 
successful  in  scaling  a  3  kW  laboratory  scale  mild 
combustion  burner  to  150  kW  level.  The  scaled 
burner  is  operated  with  two  different  fuels  and 
shown  to  achieve  mild  combustion  at  high  release 
rates  (~5.6  MW/m3)  with  both  air  and  fuel  at 
ambient  temperature.  The  design  of  its  features 
has  been  achieved  and  optimized  through  preli¬ 
minary  computations,  which  helped  in  revealing 
the  effect  of  secondary  air  position  and  injection 
velocity  on  the  recirculation  rate.  Recirculation 
rate  is  enhanced  from  153%  to  230%  by  appropri¬ 
ately  positioning  the  secondary  air  injection.  The 
distribution  of  temperature  and  02  mass  fraction 
in  the  combustion  chamber  is  essential  in  the  mild 
combustion  regime.  The  cumulative  behavior  of 
temperature-volume  and  02-volume  distributions 
show  that  02  mass  fraction  is  <0.15  in  96%  of  the 
total  volume,  and  temperature  is  >1000  K  in  93% 
of  the  total  volume.  The  presence  of  low  02  mass 
fraction  and  high  temperature  zone  in  most  of  the 
combustion  chamber  volume  is  indicative  of  mild 
combustion  at  larger  power  levels  with  current  ap¬ 
proach.  Scaling  of  air  and  fuel  jet  combination 
based  on  convective  timescales  (maintaining  be¬ 
low  80  ps)  is  an  interesting  observation.  It  needs 
to  be  explored  further  based  on  simple  jet  flame 
experiments  in  the  similar  conditions  that  exist 
in  mild  combustion  burners. 

The  experiments  with  LPG  and  producer  gas 
recorded  exhaust  emissions  of  NO  below  26  and 
3  ppm,  respectively.  The  CO  exhaust  emissions 
observed  are  1%  and  0.221%  with  heat  release 
rates  ^5.6  MW/m  .  The  measured  temperature 
and  02  gradients  in  radial  direction  are  moder¬ 
ately  small,  which  implies  that  combustion  is 
taking  place  in  mild  combustion  regime.  The  out¬ 
standing  low  chemical  emissions,  high  heat  release 


Comments 

Josette  Bellan,  Jet  Propulsion  Laboratory,  USA.  As 
you  mentioned,  Gutmark  and  Grinstein  already  per¬ 
formed  some  work  on  this  topic.  There  are  also  studies 


rates,  operation  with  two  fuels,  and  low  acoustic 
emission  features  strongly  indicate  the  potential 
for  successful  scaling  to  large  power  levels  and 
use  in  industrial  furnaces. 


References 


[1]  D.B.  Spalding,  Proc.  Combust.  Inst.  9  (1962) 
833-843. 

[2]  J.M.  Beer,  N.A.  Chigier,  Combustion  Aerodynamics. 
Wiley  Press,  New  York,  1972  (Chapter  7). 

[3]  R.  Weber,  Proc.  Combust  Inst.  26  (1996) 

3343-3354. 

[4]  J.P.  Smart,  D.J.  Morgan,  P.A.  Roberts,  Proc. 
Combust.  Inst.  24  (1992)  1365-1372. 

[5]  J.P.  Smart,  D.J.  Morgan,  Combust.  Sci.  Technol. 
100  (1994)  331-343. 

[6]  R.  Weber,  F.  Breussin,  Proc.  Combust.  Inst.  27 
(1998)  2957-2964. 

[7]  T.C.A.  Hsieh,  W.J.A.  Dahm,  J.F.  Driscoll,  Com¬ 
bust.  Flame.  114  (1998)  54-80. 

[8]  J.A.  Cole,  T.P.  Parr,  N.C.  Widmer,  K.J.  Wilson, 
K.C.  Schadow,  W.M.R.  Seeker,  Proc.  Combust. 
Inst.  28  (2000)  1297-1304. 

[9]  M.  Sadakata,  Y.  Hirose,  Fuel  73  (8)  (1994) 
1338-1342. 

[10]  A.D.  Al-Fawaz,  J.T.  Dearden,  M.  Hedley,  J.T. 
Missaghi,  M.  Pourkashanian,  A.  Williams,  L.T. 
Yap,  Proc.  Combust.  Inst.  25  (1994)  1027-1034. 

[11]  J.A.  Wunning,  J.G.  Wunning,  Prog.  Energy  Com¬ 
bust.  Sci.  23  (12)  (1997)  81-94. 

[12]  T.  Plessing,  N.  Peters,  J.G.  Wunning,  Proc.  Com¬ 
bust.  Inst.  27  (1998)  3197-3204. 

[13]  B.  Ozdemir,  N.  Peters,  Exper.  Fluids  30  (2001) 
683-695. 

[14]  P.J.  Coelho,  N.  Peters,  Combust.  Flame  123  (2001) 
503-518. 

[15]  K.  Sudarshan,  P.J.  Paul,  H.S.  Mukunda,  Proc. 
Combust.  Inst.  29  (2002)  1131-1137. 

[16]  M.  de  Joannon,  G.  Langella,  F.  Beretta,  A.  Cav- 
aliere,  C.  Noviello,  Proc.  Mediterr.  Combust.  Sym- 
pos.  3  (1999)  347-360. 

[17]  A.  Milani,  A.  Saponaro,  I  FRF  Combust.  J.  Article 
200101  (2001). 

[18]  R.  Weber,  S.  Orsino,  L.  Nicolas,  A.  Verlaan,  Proc. 
Combust.  Inst.  28  (2000)  1315-1321. 

[19]  M.  Mancini,  R.  Weber,  U.  Bollettini,  Proc.  Com¬ 
bust.  Inst.  29  (2002)  1155-1164. 

[20]  S.  Lille,  T.  Dobski,  W.  Blasiak,  J.  Propuls.  Power 
16  (4)  (2000)  595-600. 

[21]  H.S.  Mukunda,  S.  Dasappa,  U.  Shrinivasa,  in:  T.B. 
Johansson,  H.  Kelly  (Eds.),  Renewable  Energy- 
Sources  for  Fuels  and  Electricity.  Island  Press, 
Washington,  DC,  1993,  pp.  699-728. 

[22]  ABETS,  Biomass  to  Energy — The  Science  and 
Technology  of  the  IISc  Bio-energy  Systems.  ABETS, 
Indian  Institute  of  Science  Bangas,  India. 


in  the  literature  on  elliptic  and  triangular  jets.  We  have 
performed  ourselves  3-D  DNS  for  single  phase 
and  two-phase  flows  with  evaporating  drops  in  circular, 
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elliptic,  square,  rectangular,  and  triangular  geometries 
[1].  The  consensus  is  that  the  presence  of  curvature  in¬ 
duces  a  very  large  enhancement  in  stream-wise  vorticity, 
and  the  flow  in  this  direction  overtaking  the  “core”  flow 
produces  the  axis  switching.  Thus,  this  is  the  physics  of 
axis  switching,  not  the  “overdoing”  stated  in  the  presen¬ 
tation.  Moreover,  elliptic,  rectangular,  and  triangular 
geometries  result  in  axis  switching,  but  not  square  geom¬ 
etries,  which  only  exhibit  a  45°  rotation  [1,2].  Please 
comment. 
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Reply.  Vortex  dynamics,  including  axis  switching,  de¬ 
pends  on  the  nonlinear  interaction  among  all  three  vor¬ 
ticity  components.  The  author  does  not  believe  that  an 
enhanced  streamwise  vorticity  alone  provides  the  physics 
of  axis  switching,  contrary  to  the  comments  made.  The 
streamwise  vorticity  would  not  tell  the  difference  be¬ 
tween  a  major  and  a  minor  axes,  nor  the  difference  be¬ 
tween  a  square  or  a  rectangular  jet.  For  example,  the 
streamwise  vortex  pattern  and  strength  would  be  the 
same,  whether  the  initial  minor  axis  is  aligned  with  x-di- 
rection  or  y-direction,  but  the  spreading  rates  in  the  two 
cases  in  a  fixed  direction  (for  example,  x-direction)  are 
completely  different. 

As  discussed  in  the  paper,  the  author  believes  that 
axis  switching  is  due  to  the  different  magnitudes  of 
the  vorticity  components  along  the  major  and  minor 
axes,  which  led  to  different  entrainment  rates  and 
spreading  rates  along  the  two  axes.  The  differences  are 
ultimately  caused  by  the  aspect  ratio  and  curvature  ef¬ 
fects,  which  are  configuration-dependent.  The  role  of 
the  streamwise  vorticity  is  not  to  cause  the  differences 
in  flow  properties  along  the  major  and  minor  axes  di¬ 
rectly,  but  to  influence  the  vorticity  components  along 
the  two  axis  directions  through  vortex  stretching.  The 
reason  why  a  square  jet  “rotates”  only  45°  is  because 


the  jet  spreading  takes  place  in  a  symmetric  fashion  in 
the  45°,  135°,  225°,  and  315°-directions.  Another  such 
“rotation”  would  return  it  to  the  initial  square  position, 
so  it  does  not  seem  to  have  undergone  a  180°  rotation 
(axis  switching).  The  author  believes  that  the  present 
papers  offer  a  more  comprehensive  explanation  of  axis 
switching  than  before. 


Michael  Delichatsios,  University  of  Ulster ,  UK.  In 
buoyant  flames  the  velocity  accelerates  approximately 
as  Z  "  so  that  the  Reynolds  number  increases.  Thus, 
more  grid  points  are  needed  for  DNS  downstream  than 
near  the  source  to  assure  the  same  numerical  accuracy. 
Have  you  estimated  whether  your  chosen  grid  spacing 
is  adequate  throughout  the  field  of  your  calculations? 

Reply.  As  a  DNS,  the  grid  resolution  was  chosen  to  be 
sufficient  throughout  the  computational  domain.  The  fi¬ 
nal  grids  used  were  a  result  of  extensive  grid-independence 
tests.  Adequacy  of  grid  resolution  was  established  by 
monitoring  relevant  quantities  during  a  simulation  and 
by  checking  key  parameters  in  the  final  results. 


G.S.  Nathan ,  University  of  Adelaide,  Australia.  Your 
model  assumed  a  uniform  initial  velocity  profile.  In  a 
real  fire,  the  initial  velocity  is  determined  by  devolatiliza¬ 
tion,  which  in  turn,  depends  on  the  heat  transfer  from 
the  flame.  It  can  be  expected,  therefore,  that  the  real  pro¬ 
file  will  be  non-uniform,  with  the  degree  of  non-unifor¬ 
mity  increasing  with  aspect  ratio.  Please  comment  on 
the  effect  of  these  issues  on  your  results. 

Reply.  The  author  agrees  that  non-uniformity  of  the 
initial  conditions  is  an  intrinsic  feature  of  any  real  fire, 
which  was  not  included  in  the  study.  However,  such 
non-uniformity  is  expected  to  have  a  secondary,  though 
important,  effect  on  the  results  presented  because  the 
fire/flame  dynamics  considered  here  are  mainly  influ¬ 
enced  by  the  source  configuration. 


